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The entanglement of multi-atom quantum states is considered. In order to cancel noise due to 
inhomogeneous light atom coupling, the concept of matched multi-atom observables is proposed. 
As a means to eliminate an important form of decoherence this idea should be of broad relevance 
for quantum information processing with atomic ensembles. The general approach is illustrated on 
the example of rotation angle measurement, and it is shown that the multi-atom states that were 
thought to be only weakly entangled can exhibit near-maximum entanglement. 
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Progress in quantum communication and information 
processing depends on practical schemes for the massive 
entanglement of quantum systems. Advances based on 
the entanglement of atomic ensembles are undergoing 
rapid development, and involve independent approaches 
toward the entanglement of continuous light-atom vari- 
ables (see e.g., Ref. pj and references therein), and dis- 
crete variables 0, 0. Recent progress suggests that in 
the short term major advances in quantum information 
processing with atomic ensembles are likely. 

What makes quantum information processing with 
atomic ensembles very attractive is the relative simplic- 
ity of the experimental interface. As opposed to the 
case of single atoms interacting with single photons, the 
use of atomic ensembles does not require the technology 
of strong coupling cavity QED, and a free space inter- 
action between light and atomic ensemble is sufficient 
for massive entanglement. In the latter case the con- 
ditions for a sufficiently strong interaction reduce to a 
high on-resonance optical density of the atomic sample, 
a condition which is relatively easy to achieve experi- 
mentally. Therefore, the field of atomic ensemble entan- 
glement shows enormous promise as a practical way to 
implement quantum information processing. 

On the other hand, a well-known difficulty of the field 
has been how to realize the required symmetric character 
of the atom-photon interaction in the presence of inhomo- 
geneities in the local atom-field coupling. Inhomogeneity 
of this kind is a potentially fatal problem for applications 
of atomic ensembles to quantum information processing 
[4]- 

This apparent problem is rooted in a widely accepted 
notion that the multi-atom states have to exhibit a 
high degree of symmetry in order to be strongly en- 
tangled. Implicit in these arguments has been the 
fact that conventional entanglement measures, such as 
the spin squeezing parameter, are themselves symmetric 
H IE II II To the best of our knowledge, in all of 
the atomic ensemble entanglement work such symmetric 



The symmetry of entanglement measures is directly re- 
lated to the fact that additive N atom observables of the 



form A = X^fcLi^fc' which are symmetric with respect 



to permutation <-*■ d; of any two single particle op- 
erators satisfying identical commutations relations, were 
thought to provide a proper connection to the experi- 
mental situation. Even for modest degrees of asymmetry 
of multi-atom quantum states, the fact that N >> 1 im- 
poses severe limitations on the degree of entanglement of 
such symmetric observables 0. 

There is, however, no a priori reason that the mul- 
tiatom observables, and consequently, entanglement mea- 
sures, should exhibit this kind of symmetry. We would 
like to argue in favor of a more operational approach, 
in which the measurement procedure itself essentially 
defines the entanglement measure. In atomic ensem- 
ble entanglment experiments the multi-atom observables 
are typically non-symmetric, because the atoms are dis- 
tributed in space and couple to the local value of the 
electric field mode 0, 0, 0, H3, HH . Instead of symmet- 
ric operators like A, a related non-symmetric observable 
A = J^fcLi 9kO-k should be considered, where the are a 
set of numbers whose deviation from unity is associated 
with the inhomogeneity. We will show that there is no 
asymmetry induced decoherence in the entanglement of 
such matched observables. 

Our conjecture is that for any ensembles-based quan- 
tum information protocol, the asymmetry of the atomic 
ensemble cancels out and near-maximum entanglement 
is observed, provided every step of the ensemble-based 
interaction is matched to each other. In a sense the clas- 
sical noise associated with the observable asymmetry, or 
equivalently the atom-field coupling inhomogeneity, can 
be explicitly accounted for and subtracted out. 

While the notion of matched observables is quite gen- 
eral, for the sake of clarity we will illustrate it on a partic- 
ular atomic ensemble scheme: the sensitivity of rotation 
angle measurement. First let us recall the general ideas 
involved. The simplest kind of rotation measurement in- 
volves preparing an input quantum state of an N parti- 
cle system, passing it through an interferometer where a 
phase rotation is encoded, and then measuring an output 
observable sensitive to the phase rotation. The action of 
the quantum interferometer (photon or massive particle) 
can quite generally be understood in terms of the ro- 
tation of a pseudo angular momentum operator F, com- 
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prised of an ensemble of N spin 1/2 systems F = J3fe=i ^ k 
[28l | . The sensitivity of the rotation measurement is ev- 
idently conditioned on the quantum fluctuations of the 
input state of F. For a coherent spin state (CSS), a 
non-entangled state, the rotation sensitivity is given by 
the so-called standard quantum limit S<f> ~ 1/^fN. If in- 
stead we prepare a spin squeezed state (SSS), a massively 
entangled state of N particles, the phase sensitivity im- 
proves and can approach the Heisenberg limit 5(f) ~ 1/N 
|28|. This example rather clearly illustrates the practi- 
cal consequences of entanglement in quantum measure- 
ments. The crucial issue, then, is how to prepare the 
entangled SSS? One proposal involves allowing the spin 
ensemble to interact for a certain time with an auxiliary 
ensemble of n spin 1/2 particles s 1 ,* = l,...,n, prior to 
entering the interferometer [T^. Two essential new el- 
ements of this scheme, in addition to the action of the 
interferometer itself, should be noted. The interaction 
between the two spin ensembles is designed to prepare 
a SSS, and for this purpose a quantum non-demolition ( 
QND ) interaction is ideal. Moreover, an external mea- 
surement on the auxiliary spin is required S = J^iLi^ 1 ) 
in order to condition the input state to the interferom- 
eter. We shall also see later that entanglement of the 
atomic spin ensemble by means of a QND interaction 
with an ancilliary ensemble is not essential. One can also 
exploit entanglement between the two ensembles when 
the atomic ensemble is in a separable state. However, we 
will first concentrate on the former case. 

The discussion hides an important practical prob- 
lem, however, the difficulty of implementing the desired 
QND interaction. It is at this point that the issue of 
non-symmetric observables arises in practice. Analysis 
of the example of rotation angle measurement demon- 
strates how random noise associated with the asymme- 
try feeds into the measurement preventing a sensitivity 
approaching the Heisenberg limit. Consider two ensem- 
bles of spin 1/2 systems, the atoms f 1 , i = 1, ...N which 
are passed through the interferometer, and the photons 
s l ,i = 1, ...n which act as the auxiliary state preparation 
device. While at this stage we use the names atoms and 
photons merely as labels, we have in mind practical quan- 
tum optical systems where the angular momenta do cor- 
respond to a collection of atoms and a single mode probe 
light field, respectively. The total angular momenta are 

represented by S = Y17=i ^ ano - ^ = EitLi f k - We as- 
sume that initially the states of both of the ensembles 
are uncorrelated CSS, in which the average values of S 
and F are directed along the x-axis. 

The interaction between the atoms and photons is de- 
signed in order to entangle the ensemble of atomic spins, 
and is given by the QND interaction 

H = hQ J2 9kS l J k z =tMS z F z (1) 
i = l,fe=l 

where f2 is a frequency, the g k are dimensionless cou- 
pling weights, and F z = X^fcLi 9kfz 1S a non-symmetric 



atomic operator unless all of the gk are equal. In a re- 
alistic experimental scenario coupling weights of varying 
magnitude arise when an interaction of this kind is cre- 
ated using an off-resonant interaction between a single 
mode light field and a collection of atoms. Then, the dis- 
tribution of coupling weights gk maps out the variation 
of mode intensity seen by the spatially distributed atoms 
k = l,...N. 

If the atom and photon spin ensembles interact for a 
time t under Eq. the atomic spins evolve accord- 
ing to /*(t) = cos(g k xS z )f^ - sm(g k xS z )f^ /*(r) = 

sH9kXS z )fx+cos{g k xS z )fy, f z { T ) = fz, where x = fir. 
The individual spins have a dispersion of frequencies as- 
sociated with the distribution of weights, while the col- 
lective spin S satisfies 

S x \ ( cos( x F z ) -MxFz) 0\ f S x \ 

S y (r) = sin( X F z ) cos{ X F z ) ] [ S y . 
S z J V l) \S Z J 

The atoms are then passed into the interferome- 
ter where a phase rotation (f> IS imposed by means 
of a rotation about the y-axis, so that fa out — 
e i<f>F y ( T ) jk^ e -i<t>F y {r) j Q = Xj y jZ , To determine the 
phase <p we could measure the z component of total 
angular momentum at the output F z ^ out — F z cos(f> — 

sin^Efci (eos(g k xSz)fx ~ sin(g k xS z )f^. In fact we 

should really like to measure F ZjOUt — F z (t) = F z , out — F z 
the change in angular momentum due to the interferom- 
eter, cancelling any superfluous noise carried in the SSS 
input. The QND interaction enables us to obtain infor- 
mation about F z , since S y (r) ss xFzS x + S y , for xAi^ = 
Xy/N/A « 1. Hence, using (S x ) C ss = N/2, it follows 
that the operator 25 a (r)/(nx) « F z (2S x /n) + 2S y /(nx), 
is linear in F z with a coefficient of unity ( the coherent 
state is an eigenstate of S x , and thus it can be replaced 
by the c-number n/2). This operator is as close to F z as 
we can get by measuring a component of the auxiliary 
ensemble. Note that the observable which appears is the 
non-symmetric partner F z rather than the ( symmetric 
) angular momentum F z , and thus classical noise associ- 
ated with the distribution of weights gk feeds into mea- 
surements of the observable F' z = F z _ out — 2S y {r)/ (nx)- 
To determine the sensitivity of rotation measurements 
we now evaluate 

^=f V « A ^) 2 >/|^|) • (3) 

It is sufficient to evaluate (F' z ) correct to first 
order in 0, i.e., {F' z ) w ~(N<j)/2)(e^ g2 = 
-W/2)^^g^for^ X 2 « land£(A. 9 ) 2 « 1, 

where £ = nx 2 /4. Here we have used properties of the 
initial coherent spin state, and assumed that the weights 
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gk are gaussian distributed with mean value (g) = 1, and 
variance (A<?) 2 = (g 2 ) — 1. 

The variance ((AF' Z ) 2 ) may be computed at <f> = 0, 
and using conservation of F z it reduces to ((Ai 7 ^) 2 ) = 
((AF z ) 2 )+4C5 2 (r))/(n 2 X 2 )-2 (<F^(t)) + c.c) /(n X ). 
The first and second members have shot noise contri- 
butions which are cancelled by the third term, as de- 
signed, leaving residual quantum fluctuations. Explic- 
itly A(S 2 y (T))/(n 2 X 2 ) « (N + l/£ + N(Ag) 2 ) /4, assum- 
ing N\ 2 << 1) where the three terms represent shot 
noise, quantum fluctuations associated with the input en- 
tangled state and classical noise due to the distribution 
of weights, respectively. Hence, cancelling the atomic 
shot noise, we find [AF' Z ) 2 = (l + N£(Ag) 2 ) /(4f ), which 
gives a phase error of 

, P «[l+f(A 9 ) 2 ] _1 /2 



(4) 

In the symmetric limit Ag = 0, (g 2 ) = 1, this re- 
duces to <50(£) = e^ 2 /(N^). It takes a minimum 
value approaching the Heisenberg limit for £ = 1, where 
5(t>min = y/e/N [13 • Note that £ = 1 implies n% 2 = 4, 
while we have also assumed Ny^n 1 ! 2 « 4 to obtain 
the approximate expression for the variance. The two 
conditions are consistent provided the number of pho- 
tons greatly exceeds the number of atoms y/n >> N. In 
order to approach the Heisenberg limit, therefore, it is 
clear that we must satisfy (Ag) 2 << ■h a criterion which 
is very difficult to achieve for a macroscopic sample of 
atoms where typically 10 6 < N < 10 9 . 

The procedure we have described is based on a mea- 
surement of F z , but the QND interaction does not corre- 
late the auxiliary spin angular momentum S z directly to 
F z , rather to its non-symmetric partner F z . Hence the 
inhomogeneity noise will dominate the phase measure- 
ment error unless very strict limits on the distribution 
variance are satisfied. We have explicitly demonstrated 
the conventional wisdom that asymmetry washes out the 
entanglement. The important question is, can the entan- 
glement be distilled out if suitable multi-particle mea- 
surements are made to cancel the classical noise? Our 
answer is in the affirmative. The crucial point is that 
since the form of the QND interaction correlates the pho- 
ton spin to F z , we require a multi-particle measurement 
which gives F z also at the output of the interferome- 
ter. The difference signal can then in principle cancel 
the coupling distribution noise. This strategy is quite 
clear from the operational point of view taken here, but 
runs counter to standard theoretical measures of entan- 
glement based on symmetric operators, such as the spin 
squeezing parameter n = \f~N54>, which are implicitly 
sensitive to distribution noise. 

To determine which multi-particle measurements will 
achieve our goal, we consider the operator F z at the in- 
terferometer output. To first order in <f> 



Now we introduce a second set of auxiliary spins = 
1, n which are coupled with the atomic spin operators 
at the interferometer output in exactly the same way as 
spin S, i.e., H int = HttJ z F z , where J = X^LiJ 1 - In 
practice this ensemble represents a second pulse of light 
in the same spatial mode of the electromagnetic field, 
and we assume that it is prepared in a CSS pointing 
in the positive x direction, so that (J x ) = n/2. Thus 
after interacting for a time r with the atoms we have 
to sufficient accuracy for Nx 2 « 1 that S v (t) = S y + 
XF Z S X , and J v (t) = J y + xF z ,outJ x = Jy + xF z J x - 
^F x {t)J x . 

Let us now consider the difference signal A((f>) = 
J v {r) — S v {t), so that 

A(cp) =J y - S y + X F Z (S x ~S X )- <j> X F x (r) j x . (6) 

As the auxiliary spins are prepared in eigenstates of S x 
and J x , respectively, the variance in A((f> = 0) reduces to 
the sum of shot noise contributions from the two auxiliary 

spin systems ((Ai(0 = 0)) > = ((AJ y ) 2 ) + ((AS y ) 2 ) = 

n/2, and the weight distribution noise has cancelled out. 
The ensemble average of the differential change in signal 
at <p = measures F x (r) directly, i.e., 



-<i(</>)> =- X (l)(F x (r)) 
J 0=0 

X nN e-CIi+^Ag) 2 ]- 1 /2 

~ 4~ 7 9 "1 3 / 2 ' 

[l + £(Aff) 21 

and the corresponding phase error is given by 



(7) 



s<t>(0 



2 (i+e(A g )' 



3/2 



N 



(8) 



F z , ou t = F z -4>F x { T ) 



(5) 



For (Ag) 2 << 1, the minimum phase error approaches 
the Heisenberg limit 54> m i n = V2e/N, where the factor of 
v2 arises from the independent shot noise contributions 
of the two auxiliary spin ensembles. It is interesting, 
to note that even when Ag — > \J (g 2 ), the phase error 
still scales as 1/N. In other words this measurement 
scheme could operate close to the Heisenberg limit of 
phase measurement accuracy under the same conditions 
where a symmetric squeezing parameter predicts little or 
no entanglement. The measurements are appropriately 
matched to the non-symmetric observable of interest, and 
the corresponding phase error yields an improved, oper- 
ational measure of entanglement. 

It is possible to avoid using the 2nd photon ensem- 
ble Y,i — 1, ...,n, altogether. Instead, the first pho- 
ton ensemble s\i = 1, ...,n can be stored after the 
first QND interaction. After the phase rotation <f>, 
the photon ensemble is coupled to the atomic spins 
with interaction strength xi — ~X (this can be ef- 
fectively achieved by flipping the z-component of ei- 
ther atom or photon spin before the interaction takes 
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place) . The S£ component of the photon ensemble serves 
as the 0-dependent observable B{4>) = Si. We find 
B(<j>) = S y + (xF z ~ xcos(0)F z + X sin(<£)F x ) S x . It fol- 
lows that {B ((/>)) = -(1(0)), and {(AB((f) = 0)) 2 ) = 

i(^AA(0 = 0)J ), since now the shot noise of only 

one photon ensemble contributes. The phase accuracy 
<5 </>(£) is improved by a factor of V2 over that given 
by Eq.(8). It is straightforward to check that in this 
situation spin squeezing parameters 0. f° r the 
atomic ensemble always exceed or equal unity at all 
times even for perfectly symmetric coupling. Thus the 
atomic state is not "spin-squeezed", or entangled. Nev- 
ertheless, Heisenberg-limited performance of the atomic 
interferometer is achieved. Therefore, we argue that 
atom-photon entanglement in itself is the primary rea- 
son behind the improved accuracy of the considered 
QND-enhanced atomic interference measurement; the 
measurement-induced atomic squeezing is secondary and 
its use may be avoided. 

Our results are of importance for the synergy of the 
fields of cavity QED [23 and of atomic ensembles pj. 
Using a cavity of finesse F would effectively increase 
the stren gth of the ensemble-light interaction by the fac- 
tor of F • It has been suggested in the context of 
atomic ensembles that a ring cavity with running wave 
cavity modes and large Gaussian-waists should be used, 
so that all atoms see the same electric dipole mode cou- 
pling intensity [l2l l26l |27| . On the contrary, here we 
put forward an argument that suggests even standing- 
wave cavity modes can be used to realize a QND atom- 
light interaction. In this case ((Ag) 2 } — ► (g) 2 , but as we 
have shown, it is still possible to approach the Heisen- 
berg limit. In present-day experiments various kinds of 
decoherence would practically limit the amount of en- 
tanglement before the Heisenberg limit is reached. With 



that in mind, we could say that standing wave probe- 
fields are just as good as running- wave ones. The exist- 
ing high finesse cavities typically used in cavity QED are 
all of standing wave type, and the cavity modes have 
a sinusoidal spatial dependence along the cavity axis 

To what extent is the matching of observables feasible 
in practice? Let us consider the most demanding case of 
the standing wave cavity interacting with an ensemble of 
cold atoms. The atomic velocity would be on the order 
of a few cm/s, thus if the time interval between the light 
pulses < 10 fjs, the two measurements would be very well 
matched. For comparison, in recent experiments that 
demonstrated atomic memory effects in atomic ensembles 
pil25| the relevant time intervals were < 1/iS. In the case 
of a running-wave probe field, where the typical length of 
intensity variations would be a fraction of a centimeter, 
time intervals as long as 0.1 s would be possible. 

Our results should have broad relevance for situations 
involving multiatom entangled states. By presenting a 
strategy to eliminate the noise due to asymmetry of 
atomic states, we have significantly simplified the path 
towards use of atomic ensembles for quantum informa- 
tion processing. Quantum metrology applications such 
as atomic clocks, gravimeters, and atomic electric dipole 
moment searches would also benefit from enhancement in 
precision that entangled atomic states can provide. Be- 
yond the field of atomic ensembles, we hope that this 
work opens a discussion of multi-particle entanglement 
outside the symmetric entangled states that have been 
considered so far. 
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